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ABSTRACT 



^ ' Conduction in pristine conjugated polymers (other than polyacetylene) is by polaron hopping 

between sites corresponding to conjugation lengths. The strong increase of current / with voltage 

V observed for both emission-limited and ohmic contacts is due in large part to mobility increase 

as increasing field makes it more possible to overcome internal barriers, such as energy differences 

§ between sites. For emission-limited contacts an additional source of nonlinear increase of / with 

rrt ^ increasing V is greater ability to escpe return to the injecting electrode due to the image force. For 

S : ohmic contacts add.tiona. nonlineadty comes from space charge effects. We arc able to fit / V 

■ for electron or hole conduction in some poly(p-phenylene vinylene), PPV, derivatives with ohmic 

O contacts for reasonable values of the parameters involved. 
O 
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2. INTRODUCTION 



> • 

o\ : 

^ I Contact injection into conducting polymers has usually been described as tunneling of electrons or 



holes into a broad conduction or valence band, respectively, in the polymer through interface barriers 



^ ■ arising from the band offset between the polymer and the metal electrodes.0'i This picture must 
c3 ■ be corrected in several respects. First, the existence of broad continuous bands requires essentially 
^ ] infinite conjugation length whereas, it is generally acknowledged, typical conjugation lengths of PPV, 
. for example, are ~ 6 or 7 monomers. For the typical conjugation length the average spacing between 
^ ■ levels is ~ several kT at room temperature.!!! A more essential objection is that the carriers are 
' much more likely to tunnel into polaron levels than into the discrete levels in the "conduction band" 
^ or "valence band" because the polaron levels are lower in energy. The polaron levels available for 
■ electron injection lie below the lowest unoccupied molecular orbital (LUMO), while those available 
for hole injection lie above the highest occupied molecular orbital(HOMO). The separation of the 
polaron level from the LUMO or HOMO depends on chain length. For PPV, calculations give the 
separation as 0.15 eV,i or 0.2 eVjl for very long segments, increasing to ~ 0.7 eV for a 3-monomer 
long segment.! Although the polaron levels do not exist unless they are occupied by an electron or 
hole, tunneling into them may be facilitated by preparation of the appropriate chain deformation 
through the large zero point fluctuations characteristic of quasi-one-dimensional materials. 
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Once an electron or hole tunnels into a polaron level near the contact, under the influence of the 
electric fleld it moves toward the other contact, hopping from one conjugation length to another. 
The variation in conjugation length and the presence of defects result in a spread in energy of the 
hopping sites (diagonal disorder). Additional disorder is provided by the random orientation and 
variations in separation of the sites (off-diagonal disorder). As has been pointed out by a number 
of authors, the situation for hopping transport in a conjugated polymer is quite similar to that in 
molecularly doped polymers. Strong evidence for the similarity comes from the fact that the fleld 
dependence of the mobility /i generally observed for molecularly doped polymers,! 

/i = /ioe , (1) 

a being a constant, has now been observed for PPV and a couple of its derivatives. This fleld 
dependence is a result of disorder. An appropriate model for treating transport in these systems 
for low injection is the disorder model pioneered by Bassler and associates.0 The distribution in 
energy of hopping sites, in this case the polaron states, is taken as a Gaussian with variance 
Comparison of disorder theory with various experimental data for PPV leads to a ~ 0.1 eV for that 
case.0 

The behavior of a contact depends critically on the location of the polaron states relative to the 
Fermi energy, Ep^ of the metal. Internal photoemission measurements, such as those of Campbell et 
a/.,0 yielding the energy required to inject an electron from a metal into a polymer, give the energy 
difference between Ep and some average state in the polaron distribution. We denote this energy by 
W . In section II we discuss l-V characteristics for a case of large W ^ which means small injection, 
using the results of a Monte Carlo simulation based on the disorder model.El In section III we carry 
out calculations for ~ 0, or ohmic contacts, flrst for the case of only one contact injecting,lil 
then for both contacts injecting. We use the classical approach of Lampert for both cases.El'0 In 
the former case we obtain good agreement with the current / vs. voltage V characteristic of poly(2- 
methoxy, 5-(2'-ethyl-hexyloxy)-l, 4 p-phenylene vinylene], MEH-PPV, samples with only electrons 
injected, from a Ca contact^ for which photoemission gives W ~ 0.0 This calculation requires 
taking into account trapping of the electrons. A calculation of / vs. V for electrons only injected, 
including trapping, was carried out for poly(dialkoxy phenylene vinylene). The steep experimental 
variation of / with V was fltted by assuming a distribution of traps exponential in energy. These 
results are not meaningful, however, because mobility was taken as constant rather than varying 
with fleld according to Eq. (|l]). In later work by some of the same authors calculations were carried 
out for poly(dialkoxy PPV) samples with hole injection only, for which trapping may be neglected, 
with fi given by Eq. (^ and the contacts assumed ohmic.i These calculation give good agreement 
with experimental / vs. V despite the fact that W for the ITO contact that was used. The 
agreement supports the assumption that this contact was nevertheless not emission-limited for the 
range of flelds in the study.! 
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3. I-V CHARACTERISTIC FOR LARGE W 

Large W implies emission-limited contacts, which in turn usually means space charge effects may 
be neglected. Current variation with field may then be obtained by the usual Monte Carlo simulation 
for transport in disordered materials. We followed the treatment of Ref. 13. Taking Ep of the metal 
as the zero of the energy scale and the metal-polymer interface as x = 0, we may write the energy 
?7 of a polaron site in the electric field: 

U{x) = W - eEx - e^/AKx (2) 

where the last term represents the image force, k being the dielectric constant, taken as 3. In the 
presence of energetic disorder U{x) gives the mean energy e{x) of a polaron site at x. The energy e is 
assumed to have a Gaussian-distribution around e given by P(e, e) = (27rcr^)^^/^ exp[— (e — £:)^/2cr^]. 
The rate of carrier tunneling into a polaron state of energy e at a; is assumed to be Vm{x)e~^^'''^ , 
where the last factor represents the probability of the carrier having thermal energy e above Ep. The 
total number of carriers tunneling into polaron levels at x per second is then Vm{x) exp[—{e{x) — 
a'^/2kT)/kT]]^ It is reasonable to assume that this represents the initial population in the first 
layer of the polymer. Thus, as expected for large W, the mean energy of the initially populated 
sites is below e, specifically by the amount a'^/kT. It is apparent from this that increase in a, i.e., 
broadening the energy distribution of transport sites by disorder, increases the current tunneling 
into the polymer. This effect was found experimentally by Vestweber et al.^ who noted that the 
LED current increased with increase in a, indicating that the increased injection dominates over the 
decreased fi that also results from increase in disorder. 

Given the rate of carriers initially tunneling into the polymer and their distribution in energy, we 
can determine by Monte Carlo simulation the yield, i.e., the fraction of these carriers that escapes 
the return to the electrode and reaches the other boundary of the sample. Even in the absence of 
disorder the image force results in the great majority of carriers returning to the electrode at low 
fields. An analytic solution for small injection (space charge neglected) in the ordered case based on 
Eq. (2),lli@ gives the yield as 0.3 % in a field of 1.25x10^ V/cm at 300 K for the parameters of Ref. 
13. The results for the Monte Carlo simulation for a = were in good agreement with the analytic 
solution.El In the disordered case the carriers have to overcome the random barriers due to disorder 
as well as the image force. The result is that at fields ~ 10^ V/cm mild disorder [a = 0.08 eV) 
results in a yield smaller by a couple of orders of magnitudeEl even though as pointed out above, 
large a makes the injection larger. In addition the yield increases more strongly with increasing 
field, at least partially the result of increasing /i with increasing field. 

To calculate current vs. field it is necessary to multiply the yield by the number of carriers entering 
the polymer per second. The result of this for the parameters of Ref. 13 and a sample length of 120 
nm is a current increase by a factor of 10^ as the applied voltage goes from 1 to 20 eV.0 For PPV 
the increase obtained should be even larger because indications are that a ~ 0.1 rather than 0.08 
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eV. 



4. I-V CHARACTERISTIC FOR OHMIC INJECTING CONTACTS 

We consider first the case where only electrons are injected. Possion's equation may then be 
written 

{K/e)dE/dx = p (3) 

where 

p= {n-fi) + {rit- fit) . (4) 

Here e is the charge on the electron, n and nt the densities of free and trapped electrons, respectively, 
and n and fit their respective average values for the sample in thermal and electrical equilibrium (no 
applied voltage). Following LampertS we simplify the equation for current density J by neglecting 
the diffusion terms. For electrons only, with Eq. (1) incorporated, the steady current density is given 
by 

J = ncyUo exp(ai?^/^)_E' . (5) 

With the simplification of neglecting the diffusion current the boundary condition at the cathode 
interface .^O.E^oB l-V charaeteristies w.th only electrons ...jected were obtained by Parker 

on MEH-PPV samples with one Ca electrode and the other electrode Nd, Mg or Ca. As shown in 
Ref. 14 we were able to fit Parker's data with the calculations described above for both trap-free 
[rit = fit = ^) and all-traps-filled {rit = Nt, the total trap density) cases.lii In the latter case there 
is a threshold for current flow because the repulsion of the charges in the traps prevents injection of 
additional charge at low fields. However, beyond the threshold, above ~ 5 V or 4x10^ V/cm for the 
particular samples, there is a good fit to the data under the assumption that all traps are filled. For 
the trap-free case the parameters for the fit were /io = 5x 10^^^ cm^/Vs, a = 8x 10^^ cm^/^/V^/^. For 
the trap-filled case the best fit was obtained with /io = 5 x 10^^ cm^/Vs, a = 4.5 x 10^'^ cra^^"^ /V^^"^ 
and Nt = 10^''/cm^. (It is reasonable to assume fi = fit = 0). The latter set of parameters, 
particularly /xq, appears to be in better agreement with values determined by other methods. Hole 
mobility vs. field, measured for PPV by observing the time delay after injection for the appearance 
of luminescence, yields by extrapolation /io = 5 x 10^^ cm^/Vs.i Electron mobility in PPV and 
its derivatives is generally found at not too high fields to be a couple of orders of magnititude 
smaller than hole mobility, /ip, the differences attributed to trapping, by carbonyls for example.0 
It is reasonable that /i„ increases as E increases and thus trap-filling increases. The finding that 
/^n ~ /^p in the trap-filled limit is not unexpected because bandstructure is similar for electrons 
and holes in PPV.0 It is interesting to note that /ig is different for different PPV derivatives. For 
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dialkoxy-PPV, for example, Blom et a/ obtained fio = 5 x 10~^ cm^/Vs for holes.0 We suggest that 
the differences reflect different degrees of disorder. Support for this is the fact that in dialkoxy-PPV 
/i can be measured by tranditional time-of-flight technique, in which the arrival time at the other 
electrode of an injected pulse of carriers is observed, whereas transport in PPV and MEH-PPV is 
too dispersive for that method to work. The value of a found by Karg et al. is 6 x 10~^ cm^/^/V^/^,0 
close to the value of 4.5 x 10"^ cm^/^/V^/^ found here, and also to the value 5.4 x 10"^ cm^/^/V^/^ 
found for dialkoxy-PPV.^ The value of Nt for the particular sample used is, of course, unknown but is 
similar to the value estimated by Campbell et al. for MEH-PPV from the change in capacitance with 
forward bias.!! From the good fit at high voltage we conclude that the current there is space-charge 
limited. It does not vary as V"^ because of the strong field-dependence of /x. 

We consider now the case of both electrons and holes injected. For this analysis we simplify 
Possion's equation by taking 

p = n-p + Nd, (6) 

where A^^^ is the density of traps filled by electrons. For the current density we use Eq. with /io 
replaced by , and a similar term added for holes with /i° replaced by /i^. We assume a is the same 
for electrons and holes. Steady state requires that 

dJ djn djr, 
dx dx dx 

where = nuJE^E and ?„ = piiJE)E. For the steady state, on the simplifying assumption that 
the lifean. fo. Lib r^l constant. WependLt of we obtain! 

djn/dx = -n{x)/Tn , (8) 
djp/dx = n{x)/Tn . (9) 

Eliminating n and p from Eqs. (|^), (|), (|]) and (|^), and using Eq. (|l]), we obtain 

Ee'^±(Ee'^-f) - -^Ee'^f = -^^^ . (10) 

dx dx fx^Tn dx efx^n^Tn 



This equation can be solved by assuming 



^Ve d d 



^g^v^ - = ^ . (11) 
dx dy 

The solution is 

Eiy) = C^e''y-iA/b)y + C2 (12) 

where b = 1/fipTn, A = —J/e/ipii^Tn and Ci, C2 arbitrary constants to be determined by the 
boundary conditions. 
We assume both contacts are ohmic, so the condition at the cathode is 
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E{y = 0) = = 0) = 



(13) 



and at the anode 

E{x = L) = E{y = ye) = (14) 

where t/c is determined from 

L = r E(y)e''^/^^dy . (15) 
Jo 

With these boundary conditions 

Ci = {A/b)y,{e'^^ - C2 = -Ci . (16) 

The apphed voltage is given by 

V= 1^ E{x)dx= r E\y)e''y^Uy . (17) 

These equations have been used to fit data for an MEH-PPV diode with one Ca contact and the 
other contact Au or ITO, all of which may reasonably be considered ohmic. Shown in Fig. 1 is the 
fit for N^ — O and the parameters for that fit. It should be noted that the //^ value is similar to that 
obtained earlier for the trap-free single carrier case. To further elucidate the meaning of the results 
we calculated jn and jp, the electric field, n, p and the ratio p/n as functions of x for a particular 
current density, J = 0.2 mA/mm^ corresponding to V — 7.4 V. It was found that the large ratio 
fj>p/ /J'n causes p » n over 2/3 of the sample. This results in jp >> j„ except in the immediate 
vicinity of the cathode. To maintain constant current given the large ratio of /ip to requires that 
the electric field be larger near the cathode than the anode, rising more rapidly from £■ = at the 
cathode than at the anode. If we allow Ndy^O and /ip — as expected for the trap-filled case, the 
curves become more symmetrical in n and p and the field is more symmetrical about the middle of 
the sample. 

5. CONCLUSIONS 

Because of the short conjugation lengths, characteristic of PPV and other conducting polymers, 
there are no broad continuous bands and carriers are most likely to tunnel into polaron levels with 
little or no hindrance by a barrier. Random variation in conjugation lengths, and thus site energies, 
as well as other sources of disorder, result in the hopping mobility of polarons increasing exponentially 
with y/E, as has been shown exponentially for several PPV derivatives. This increase of // with E 
contributes to the strong increase of / with V observed for LEDs, for example. We have studied the 
variation of / with V for polymer samples with emission-limited contacts and with ohmic contacts. 
In the former case we did not assume any /i variation with E but relied on the usual Monte Carlo 
simulation to determine first the number of carriers that make it across the sample to the other 
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electrode as a function of E, and then the current. We find that / increases steeply as V increases, 
both because carriers can better overcome the image force and more easily surmount barriers due 
to random site energy variations. 

In the case of ohmic contacts, space charge effects also contribute to the steepness of / vs. V. It is 
now documented that the / vs. V characteristic is well fitted by the use of Poisson's Equation and 
the current continuity condition, with insertion of the /i dependence on E, for hole only transport in 
dialkoxy PPV.I^ Fits to Parker's / vs. V data for hole transport only and electron transport only in 
MEH-PPV,0 suggest that trapping is important in both those cases. A good fit has been obtained for 
high enough fields that all traps should be filled, with parameters that are in reasonable agreement 
with the values found from other experiments. More work is needed to fit / vs. V for the case where 
both carriers contribute. 
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Fig.l. Fit of theory (solid curve) to data (dots) of Ref. 1 (Fig. 13) for MEH-PPV with electrons 
injected at Ca contact, holes injected from ITO. The parameters for the fit are = 0, a — 7 x 10~^ 
cmVVVi/2, = 10-9 cmVVs, /x^ = 0.01/^°, r - 1.4 x 10"^ s. 
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